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Abstract 

A theorem of Bayer and Teichmann [BT] implies that if a finite real multisequence 
/3 = /J^^'^) has a representing measure, then the associated moment matrix Md admits 
positive, recursively generated moment matrix extensions Md+i, Md+2, ■ ■ ■■ For a bi- 
variate recursively determinate Md, we show that the existence of positive, recursively 
generated extensions Md+i, ■ ■ ■ ,M2d-i is sufficient for a measure. Examples illustrate 
that all of these extensions may be required to show that /3 has a measure. We describe 
in detail a constructive procedure for determining whether such extensions exist. Under 
mild additional hypotheses, we show that Md admits an extension Md+i which has many 
of the properties of a positive, recursively generated extension. 
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1 Introduction 

Let /3 = /3(2'^) := {ft^ji^ 

j>o,i+j<2d denote a real bivariate moment sequence of degree 2d. 
The Truncated Moment Problem seeks conditions on /3 for the existence of a positive Borel 
measure on such that 



/ xVd/i (i,j>0, i + j<2d). (1.1) 

7IR2 



A result of jCF5j shows that /3 admits a finitely atomic representing measure fi (as in (jl.ip ) 
if and only if = Md{/3), the moment matrix associated with /3, admits a flat exten- 
sion Md+k+i, i-e., an extension to a positive semidefinite moment matrix Md+k+i such that 
rank M^+k+i = rank M^+k- The extension of this result to general representing measures 
follows from a theorem of C. Bayer and J. Teichmann |BTj . which implies that if /3 has a 
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representing measure, then it has a finitely atomic representing measure (cf. |FH Section 2], 
[CFMl Section 1]). At present, for a general moment matrix, there is no known concrete test 
for the existence of a flat extension M^+k+i- In this note, for the class of bivariate recursively 
determinate moment matrices, we present a detailed analysis of an algorithm of |Flj that can 
be used in numerical examples to determine the existence or nonexistence of flat extensions 
(and representing measures). This algorithm determines the existence or nonexistence of 
positive, recursively generated extensions M^^i, . . . ,M2d-i, at least one of which must be 
a flat extension in the case when there is a measure. Theorem 12.51 shows that there are 
sequences (3^'^^^ for which the first flat extension occurs at M2d-i, so all of the above exten- 
sions must be computed in order to recognize that there is a measure. This result stands in 
sharp contrast to traditional truncated moment theorems (concerning representing measures 
supported in M, [a, b], [0, +oo), or in a planar curve of degree 2), which express the existence 
of a measure in terms of tests closely related to the original moment data (cf. Remark 12.61 
below and [CFl] . [CF2] . [CF3] . |CFMj . [FT]). Here we see that, at least within the frame- 
work of moment matrix extensions, we may need to go far from the original data to resolve 
the existence of a measure. In Theorems 12.31 and 12.131 we show that under mild additional 
hypotheses on Md, the implementation of each extension step, from Md+j to Md+j+i, leading 
to a flat extension Md+k+i^ consists of simply verifying a matrix positivity condition. 

Let Vd = M[a;, y]d denote the bivariate real polynomials of degree at most d. For p € Vd, 
p{x,y) = aijx'^y^ , let p := (uij) denote the vector of coefficients with respect to the 

i,j>0,i+j<d. 

basis for Vd consisting of monomials in degree-lexicographic order, i.e., 1, x, y, x'^,xy, y^, . . . , x'^, 
. . . ,y'^. Let Lp : V2d — > ^ denote the Riesz functional, defined by Lis{ aijX^y^) := 

''^^^aijPij. The moment matrix Md, whose rows and columns are indexed by the monomials 

in Vd, is defined by {MdP,q) := Li3{pq) {p,q G Vd)- We denote the successive rows and 
columns of Md by 1, X, y, . . . , X'^, . . . , thus, the entry in row X^Y^ , column X^Y^, 
which we denote by {X''Y^,X'Y^), is equal to (3i+k,j+i- We may denote a linear combi- 
nation of rows or columns by p{X, Y) := ^ aijX^Y^ for some p = ^ aijX^y^ G Vd', note 
that p{X, Y) = Mdp. We say that Aid is recursively generated if ker Md has the following 
ideal-like property: 

P,q,pq^Vd, p(x,y) = o^(pg)(x,y) = o. (1.2) 

If /? has a representing measure, then Md is positive semidefinite and recursively generated 
|CF5j (and in one variable these conditions are sufficient for the existence of a representing 
measure |CF1] ). Moreover, from |BTj . /3 actually admits a finitely atomic representing mea- 
sure which therefore has finite moments of all orders; it follows that Md admits positive, re- 
cursively generated moment matrix extensions of all orders, namely Md+iln], ■ ■ ■ , Md+k[lj\i ■ ■ ■■ 
Let us consider a moment matrix extension 

M, B{d+l) 
^''+'=[Bid+lf C{d+1) 

where the block B{d + 1) includes new moments of degree 2d -\- 1 (as well as old moments 
of degrees d + 1, . . . ,2d), and block C{d + 1) consists of new moments of degree 2d + 2. 
We denote the columns of B{d + 1) by . . .,Y'^+^, and we say that {Md B{d + 1)) is 
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recursively generated if (jl.2|) holds in its column space, but with p,q,pq £ Vd+i- Afrf+i is 
positive semidefinite if and only if (i) is positive semidefinite; (ii) Ran B{d+1) C Ran 
(equivalently, B{d + 1) = MdW for some matrix W); (iii) C{d + I) h := W^MdW (cf. 
[CF2j ). (Here and in the sequel, for a real symmetric matrix A, we will write A ^ (resp. 
A y 0) to denote that A is positive semidefinite (resp. positive semidefinite and invertible).) 
If Md+i ^ 0, then we also have (iv) each dependence relation in Col (the column space of 
Md) extends to Col M^+i. In the sequel we say that M^+i is an RG extension if properties 
(i), (ii), and (iv) hold and M^+i is recursively generated (so, in particular, {Ma B{d+ 1)) 
is recursively generated). In the sequel, we provide sufficient conditions for RG extensions; 
note that to verify that an RG extension is positive semidefinite and recursively generated, 
it is only necessary to verify condition (iii). 

For a general M^, a significant difficulty in determining the existence of a flat extension 
Md+k+i is that there may be infinitely many positive and recursively generated extensions 
Md+i. If one such extension does not admit a subsequent flat extension, this does not 
preclude the possibility that some other extension does. In the sequel, we focus on the class 
of recursively determinate moment matrices (RD) introduced in |F1| (cf. |F2j ). These are 
characterized by the property that there can be at most one positive, recursively generated 
extension, and there is a concrete procedure (described below) for determining the existence 
or nonexistence of this extension. Since such an extension is also recursively determinate, we 
may proceed iteratively to determine the existence or nonexistence of positive and recursively 
generated extensions 

Mrf+i,...,M2d-i. (1.3) 

As we discuss below, the existence of the extensions in (jl.3p is equivalent to the existence of 
a flat extension M^+fc+i and, in fact, one of the extensions in (jl.3p is a flat extension of M^. 
(If Mj^i is positive semidefinite, then Mj is positive semidefinite and recursively generated 
|CF5j . so, using also jBTj . it follows that ()1.3p is equivalent to the existence of a positive 
semidefinite extension M2d-) 

A bivariate moment matrix Md admits a block decomposition Md = (i?[i, j])o<.tj<d) where 



B[i,3] 



\ ••• /3o,i+j/ 

Thus, B[i,j] is constant on each cross-diagonal; we refer to this as the Hankel property. Note 
that in the extension Md+i^ B{d + 1) = {B[i^d + l])o<i<d, and all of the new moments of 
degree 2d+\ appear within block B[d,d+ 1], either in column X'^'^^ (the leftmost column) 
or in column (on the right). Similarly, all new moments of degree 2d + 2 appear in 

column X'^'^^ or column y^+i of C{d + 1) (= B[d + l,d + 1]). In the sequel, by a column 
dependence relation we mean a linear dependence relation of the form X^Y^ = r(X, 1"), where 
deg r < i + j and each monomial term in r strictly precedes x^y^ in the degree-lexicographic 
order; we say that such a relation is degree reducing if deg r < i + j. A bivariate moment 
matrix Md is recursively determinate if there are column dependence relations of the form 

X"" = p{X,Y) {p G Vn-i. n<d) (1.4) 

and 

Y"^ = q{X, Y) {q G q has no term, m < d), (1.5) 
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or with similar relations with the roles of p and q reversed. In the sequel, we state the main 
results (Theorems 12.31 and l2.5p with p and q as in (ll.4p - ()1.5p . but these results are valid as well 
with the roles of p and q reversed. In Section[2]we show that if is recursively determinate, 
then the only possible positive, recursively generated (or merely RG) extension is completely 
determined by column relations X'^+^ = (a;'^+i-»(X, and V^^^ = {y'^+'^-"^q){X,Y). 

The most important case of recursive determinacy occurs when is positive and flat, 
i.e., rank = rank Mfi-i (equivalently, each column of degree d can be expressed as a linear 
combination of columns of strictly lower degree) . A fundamental result of |CF2| shows that 
in this case admits a unique flat extension M^j^i (and a corresponding rank M^-atomic 
representing measure). In this paper, we stay within the framework of recursive determinacy, 
but relax the flatness condition, and study the extent to which positive, recursively generated 
extensions exist. 

Our main results are Theorems 12.31 and 12.131 which give sufficient conditions for RG 
extensions, and Theorem 12.51 which shows that the number of extension steps leading to a 
flat extension is sometimes proportional to the degree of the moment problem. Theorem 
12.31 shows that if is positive and recursively generated, and if all column dependence 
relations arise from (|1.4p or (11. Sp via recursiveness and linearity, then admits a unique 
RG extension. In general, this extension need not be positive semidefinite (see the discussion 
preceding Example II. ip . but \i d = n + m — 2, then this extension is actually a flat extension, 
so (3 admits a representing measure (Corollarv l2.2p . Additionally, we show in Theorem 12. 131 
that if Md is positive semidefinite, recursively generated, and recursively determinate, and 
if all column dependence relations are degree-reducing, then again admits a unique RG 
extension. However, we show in Example 1 2 . 1 2 1 that if all of the column relations are degree- 
reducing except that deg q = m, then need not even admit a block B{d + 1) consistent 
with recursiveness for (M^ B{d + 1)). In Theorem 12.51 we show that for each d, there exists 
/3 = with Md(/3) e RD, such that in the sequence of positive, recursively generated 

extensions, M^+i, . . . ,M2d-i, the first flat extension is M2d-i, so the determination that a 
measure exists takes the maximum possible number of extension steps. Moreover, at each 
extension step, M^+j satisfies the hypotheses of Theorem 1 2. 3 ^ so it is guaranteed in advance 
that the next extension M^+i+i is well-defined and recursively generated; only its positivity 
needs to be verified. In general, however, the existence of a positive, recursively generated 
extension M^+i does not imply the existence of a measure. In Section [3] we answer \F1\ 
Question 4.19] by showing that if, under the hypotheses of Theorem 12.3^ Mj, does admit a 
positive, recursively generated extension M^+i, then M^+i may also satisfy the conditions of 
Theorem 12.31 but need not admit a positive, recursively generated extension Md+2, and thus 
Md may fail to have a measure. 

We conclude this section by reviewing and illustrating [Fll Algorithm 4.10] concerning 
extensions of recursively determinate bivariate moment matrices. We may assume that Md 
is positive and recursively generated, for otherwise there is no representing measure. (We 
note that in numerical problems, positivity and recursiveness can easily be verified using 
elementary linear algebra.) To define block B{d + 1) for an extension Md+i, note that 
blocks B[0,d + 1], . . . , B[d — l,d + 1] consist of old moments from Md- To define moments 
of degree 2d + 1 for block B[d,d + 1], we first use (II. 4p and recursiveness to define the "left 
band" of columns, x'^+V'^+i-^-" := (a;'y'^+^-'~»(X, F) {0 < i < d + 1 - i). In block 
B[d,d + 1], certain "new moments" in column X"'Y'^^^^"' can be moved up and to the right 
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along cross-diagonals until they reach row X'^ (the top row of B[d,d+ 1]) in columns of the 
"central band," x"-iy'^+2-", . . . , r^^^g^ ^^^^^^ ^^^^ ^^^^ ^^^^^ 

l^-^d+i-mYm^X'^) (the entry in row X'^, column X'^+^-^F'") by means of 

l^^d+l-mym^ X'^) := ( (x'^+^-'^q) (X, Y) , X'^) . (1.6) 

This value may be moved one position down and to the left along its cross-diagonal and then 
used to define I^X'^+^-'^y'^ , X'^-^Y) := ((x'^+i^™g)(X, y), X''-iy). We repeat this process 
successively to complete the definition of column X'^~^^~"^Y"^ in B[d,d + 1] as well as the 
definition of the central band of columns in this block. We next complete the definition of 
B[d,d+ 1] by successively defining the "right band" of columns, X'^~'^Y"^^^ , . . . , Y'^^^, using 

^d+l-m-iym+i .- (x'^+^-''^-iyiq)(X , Y) {0 < i < d + 1 - m) . 

It is necessary to check that the values in the central and right bands, as just defined, 
are compatible with values in the left band, and, more generally, to verify that B(d + 1) is 
a well-defined moment matrix block. If this fails to be the case, there is no measure. If 
B{d + 1) is well-defined, we next check that Ran B[d -|- 1) ^ Ran M^, for if this is not the 
case, then there is no measure. Assuming the range condition is satisfied, ()1.4p and ()1.5p 
will hold in the columns of B{d + 1)"^ (the transpose). We then apply recursiveness and the 
method used just above in defining B\d,d+\\ to attempt to define C{d+ 1) = B[d+ l,d+ 1]. 
Assuming that C{d+1) is well-defined, we further check that M^+i is positive and recursively 
generated. If any of the preceding steps fails, there is no representing measure. Our main 
results (Theorems 12.31 and I2.13P show that if all column relations come from (jl.4p or (jl.Sp via 
recursiveness and linearity, or if (|1.4p - (jl.Sp hold and all column dependence relations are 
degree-reducing, then all of the preceding steps are guaranteed to succeed, except possibly 
the positivity of M^+i; thus M^j^i is at least an RG extension. 

If Mfi+i, as just defined, is positive and recursively generated, then, since it is also recur- 
sively determinate, we may apply the above procedure successively, in attempting to define 
positive and recursively generated extensions M(i+2 ) ^d+'i ) • • • • Note that the central band of 
degree d in has n+m—d—1 columns, x^^^y^^""*"^, . . . , x'^'^^~'^Y"^~^ . In each successive 
extension M^j^k-, the number of columns in the central band of degree d+k is n + m — d—l — k. 
Thus, after at most n + m — d — 1 extension steps, either the extension process fails, and 
there is no measure, or the central band disappears and there is a flat extension, at or before 
M„_i_m-i, and a measure. (Note that since n,m < d, this refines our earlier assertion that 
a fiat extension occurs at or before M2d-i-) Another estimate for the number of extension 
steps is based on the variety of M^, defined as V = V{Md) ■= 2r, where Zr is 

the set of real zeros of r. It follows from |F1] that the number of extension steps leading to a 
fiat extension is at most 1 -|- card V — rank M^. Note also that when a measure exists, it is 
supported inside V |CF5j . so its support is a subset of the finite real variety determined by 
— p(x, y) and — q{x, y). 

Examples are known where the RG extension M^+i is not positive semidefinite (cf. |Flt 
Example 4.18], [CFM, Theorem 5.2], both with n = m = d = 3, and the example of Section [3] 
(below), with d = 5, n = m = 4). We next present an example, adapted from |F2^ Example 
5.2], which illustrates the algorithm in a case leading to a measure. 
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Example 1.1. Let d = 3 and consider 
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We have M3 ^ 0, M2 ^ 0, and rank M3 
rank M3 = 8, then the two column relations 
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8 d = 2026881 - 2844c + c^. When 

e 



and 

Y'' = q{X,Y), 

where q{x, y) := (5715 - 4c)x + 10(-1428 + c)y - 3(-2853 + 2c)x^y + (-1422 + c)xy^. Let 
ri(x,y) = y — x^ and r2{x,y) = y^ — q{x,y). With these two column relations in hand, 
Theorem 12.31 guarantees the existence of a unique RG extension M4. To test the positivity 
of M4, we calculate the determinant of the 9x9 matrix consisting of the rows and columns of 
M4 indexed by the monomials 1, x, y, x'^,xy, y^, x^y, xy'^,x'^y'^. A straightforward calculation 
using Mathematica shows that three cases arise: 

(i) c < 1429: here M4 0, so M3 admits no representing measure; 

(ii) c = 1429: here M4 is a flat extension of M3, so by the main result in [CF2j . M3 admits 
an 8-atomic representing measure; 

(iii) c > 1429: here M4 is a positive RG extension of M3 with rank 9. Although M4 is not 
a flat extension of M3, it nevertheless satisfies the hypotheses of Theorem l2.3t so Corollary [23] 
implies that M4 admits a flat extension M5, and therefore M3 has a 9-atomic representing 
measure. Moreover, since the original algebraic variety V = ViM^) associated with M3, 
Zrj^C\2r2, can have at most 9 points (by Bezout's Theorem), it follows that V = V{M^). 
This algebraic variety must have exactly 9 points, and thus constitutes the support of the 
unique representing measure for M3. 

To illustrate this case, we take the special value c = 1430, so that q{x,y) = —5x + 
20y — 21x^y + Sxy"^. Let a := ^\/o — 2a/5 and 7 := \/5q. A calculation shows that 
V = {{xi,xl)}^^^, where xi = 0, X2 = ^(-1 - V5) ^ -1.618, X3 = ^(1 - Vs) w -0.618, 
X4 = — X3 0.618, x^ = —X2 ~ 1.618, xq = —a — 7 —1.176, xy = —a + 7 w 0.449, 
xg = —xj w —0.449 and xg = — xg ~ 1.176. M3 satisfies the hypothesis of Theorem 12.31 with 
n = m = 3, so we proceed to generate the RG extension M4. This extension is uniquely 
determined by imposing the column relations = XY, X^Y = Y^, XY^ = {xq){X, Y), and 
y4 = (yq)(X,Y) (first in ( M3 B{A) ), then in ( B{Af C(4) )). A calculation shows 
that, as expected, these relations unambiguously define a positive moment matrix M4 with 
rank M4 = 9 (> 8 = rank M3). It follows that M3 admits no fiat extension M4, so we proceed 
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to construct the RG extension M5, uniquely determined by imposing the relations = X^Y , 
X^Y = XY^, X^Y^ = Y^, X2y3 = lx^q){X,Y), XY^ = {xyq){X,Y), Y^ = {y^q){X,Y). 
A calculation of these columns (first in ( M4 -6(5) ), then in ( -6(5)'^ ^"(5) )), shows 
that, as again expected, they do fit together to unambiguously define a moment matrix M5. 
From the form of q{x,y), we see that M5 is actually a flat extension of M4, in keeping 
with the above discussion. Corresponding to this flat extension is the unique, 9-atomic, 
representing measure // = fiM^ as described in |CF5j . Clearly, supp fi = V, so n is of the 
form /i = Yll=i Pi^{xi,x^)- To compute the densities, we use the method of |CF5j and find 
Pi = l = 0.2, P2 = P5 = « 0.069, p3 = p4 = « 0.181, P6 = P9 = « 0.131, 

and pr = ps = ^^q^t^ ~ 0.019. Thus, the existence of a representing measure for /S*-^^ is 
established on the basis of the extensions M4 and M5, in keeping with Theorem l2.31 Note that 
in this case, the actual number of extensions leading to a flat extension can be computed as 
either n + m — d—1 or as l + card V — rank M3, which is consistent with our earlier discussion. 

□ 

Acknowledgment. Examples 11.11 and 12.121 and the example of Section [3l were obtained 
using calculations with the software tool Mathematica fW^ . 



2 The extension of a bivariate RD positive moment matrix 

In Theorem 12.31 (below) we show that a positive recursively determinate moment matrix 
Mrf = Mrf(/3), each of whose column dependence relations is recursively generated by a 
relation of the form 

X^ =p{X,Y), (peVn-i) (2.1) 

or 

Y"' = q{X,Y), (qeVm), (2.2) 

(where n,m < d are fixed and q has terms with v < m), always admits a unique RG 
extension 

Md B{d+1)\ 
,B{d+lf C{d+l)J- 

The main step towards Theorem 12.31 is the following result, which shows that (as above) 
admits an extension block B(d + 1) that is consistent with the structure of a positive, recur- 
sively generated moment matrix extension M^+i. 

Theorem 2.1. Suppose the bivariate moment matrix Mrf(/3) is positive and recursively gen- 
erated, with column dependence relations generated entirely by h2. 1\) and h2. via recursive- 
ness and linearity. Then there exists a unique moment matrix block B{d + 1) such that 
(Mrf B{d + 1)) is recursively generated and Ran B{d + 1) C Ran M^- 

The hypothesis implies that the column dependence relations in Md are precisely those 
of the form 

X^+iyi = {x'y^p) (X, y ) (i , J > 0, i + J + n < d) (2.3) 

and 

X^y™+' = (xVg)(^,>") {k,l>Q,k + l + m<d). (2.4) 
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In particular, the degree d columns X'^, . . . ,X"'Y'^ " are recursively determined in terms of 
columns of strictly lower degree. Since, by dMD, each column x'^-'"-'=y"^+'= {0 < k < 
d — m) may be expressed as a linear combination of columns to its left, it follows that if 
n < d — m + 1, then is flat. Since a flat positive moment matrix admits a unique positive, 
recursively generated extension (cf. |CF2j ). we may assume that not every column of degree 
d is recursively determined, i.e., n > d — m + 1, or 

n + m> d+l. (2.5) 

We may denote 

X^=p{X,Y)^ arsX^Y' (2.6) 

r,s>0,r+s<n—l 

and 

Y^ = qiX,Y)^ Yl KvX^Y^ 

u,v>0,u+v<m,v<m 

Y aabX''Y''+ Y lceX'Y\ (2.7) 

a,b>0,a+b<m—l c,e>0,c+e=m,e<m 

Thus, in any positive and recursively generated (or merely RG) extension M^^+i, certain 
columns of B{d + 1) are recursively determined. On the left of B{d + 1), there is a band of 
columns, 

= Y arsX^^+^Y'+'^+^-^'-f (0 < / < d + 1 - n) (2.8) 

r,s>0,r+s<n—l 

each of which is well-defined as a linear combination of columns of M^^. On the right of 
B{d + 1) there is another band of recursively determined columns, 

^d+l-m-gym+g ._ ^^d+l-m-g yg ^^^^ ^y) 

= Y 6™X"+'^+i-"^-5y"+f (o<g<d+l-m). (2.9) 

u,v>0,u+v<m,v<m 

If deg q = m, the sum in (12. 9p may involve columns from the middle band, x^~^'^~^^~'^~^Y^~^^ 
{u + V = m,u + d + 1 — m — n < g < m — v), which has not yet been defined, so some care is 
needed in implementing (12. 9p . 

The proof of Theorem 12.11 entails two main steps, which we prove in detail in Section 
the construction of the block B{d + 1), and the verification of the inclusion Ran B{d + 1) C 
Ran Mfi- Assuming that we have already built a unique block B[d + 1) consistent with the 
existence of a positive, recursively generated extension 

^^^^+^^\B{d+lY C{d + l] 

we next use this to construct a unique block C{d+1) consistent with the existence of an RG 
extension. 
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Corollary 2.2. // satisfies the hypotheses of Theorem \2.1l then there exists a unique 
moment matrix block C = C(n + 1) consistent with the structure of an RG extension Md+i- 

Proof. In any RG extension M^+i the column relations (|2.8p and (j2.9p must hold. The proof 
of Theorem 12.11 shows that these relations define a unique moment matrix block B = B[d-\- 1) 
consistent with positivity and recursiveness. To define C = C(n + 1), we may formally repeat 
the proof of Theorem 12.11 concerning the well-definedness and uniqueness of block B{d+ 1), 
but applying the argument with replaced with B{d + 1)^, and B[d, d + 1] replaced with 
G = B[d + 1, d + 1]. In brief, we use B{d + 1)^ and (12. 8p to define the left recursive band in 
G. We then define column X'^^^~"^Y"^ by applying (j2.9p successively, starting in row 
so that this column is Hankel with respect to the central band, which we are completing 
simultaneously. We then use (|2.9p to successively define the remaining columns on the right. 
Lemma 14.31 can be used to show that the left band is internally Hankel, and an adaptation 
of the argument in Lemma 14.41 can be used to show that column X'^~^^~'^Y"^ is Hankel with 
respect to the left and central blocks. Finally, the argument of Lemma 14.51 can be adapted 
to show that that the right band is also Hankel. □ 

By combining Theorem 12.11 with Corollary 12. 2^ we immediately obtain the first of our 
main results, which follows. 

Theorem 2.3. If is positive, with column relations generated entirely by 12. 1\} and \2. 2\) 
via recursiveness and linearity, then admits a unique RG extension M^+i, i.e.. Ran B[n+ 
1) C Ran Md, i2.8\ )- [2lV\} hold in Gol M^^i, and M^+i is recursively generated. 

Corollary 2.4. If satisfies the hypotheses of Theorem \2.3\ and d = n + m — 2, then 
admits a flat moment matrix extension M^^i (and /3 admits a rank M^-atomic representing 
measure). 

Proof. Each column in the left band is, from (12. Sp . a linear combination of columns of strictly 
lower degree. Since d = n + m — 2, there is no central band in the construction oi B{d+ 1) in 
Theorem 12.11 and of C((i+ 1) in Corollary 12. 2[ It thus follows from (j2.9p that each column 
in the right band is also a linear combination of columns of strictly lower degree, so M^+i is 
a flat extension. □ 

To illustrate Corollary 12.41 in the simplest case, let n = m = d = 2 and suppose that M2 
satisfies the hypotheses of Theorem 12. 31 It follows from |CF4j that M2 admits a representing 
measure if and only if the equations x'^ — p{x, y) = and — q{x, y) = have at least 4 
common real zeros. Corollary 12.41 implies that the latter "variety condition" is superfluous; 
indeed, from Corollarv 12.41 there is a representing measure, so [CF4j implies that the system 
must have at least 4 (= rank M2) common real zeros. 

Note that if Mfi{(3) satisfies the hypothesis of Theorem 12.31 then the existence or nonex- 
istence of a representing measure for /? will be established in at most d — \ extension steps 
(after which the central band would vanish and every column of M2d~i would be recursively 
determined). The next result shows that for every d > 2, there exists Md{/3), satisfying 
the conditions of Theorem 12.31 for which the determination that a representing measure ex- 
ists entails the maximum number of extension steps, each of which falls within the scope of 
Theorem ESI 
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Theorem 2.5. For d > \, there exists a moment matrix Md, satisfying the conditions of 
Theorem \2.3l for which the extension algorithm determines successive positive, recursively 
generated extensions M^+i, ■ ■ ■ ,M2d-i, and for which the first flat extension occurs at M2d-i- 
Moreover, each extension M^+i satisfies the conditions of Theorem \2.3\ so to continue the 
sequence it is only necessary to verify that the RG extension M^_|_j_|_i is positive semidefinite. 

Remark 2.6. To illustrate the significance of Theorem [231 l^t us compare it to the following 
result of |CF6[ Theorem 1.2]: If Md{(3) is a bivariate moment matrix with a column relation 
p{X,Y) = {deg p < 2), then /3 has a measure if and only if Md is positive, recursively 
generated, and rank Md < card V{Md). In this result, we see that the existence of a measure 
can be determined directly from the data by establishing the positivity, rank, and variety 
of Md- By contrast, in Theorem 12.51 we see that it may be necessary to extend Md to 
M2d-i in order to establish that a measure exists. In this sense, within the framework of 
moment matrices, we see that the general case of the truncated moment problem cannot be 
solved in "closed form." We may therefore seek to go beyond the framework of moment 
matrices. Recall that for /3 = is positive if p G V2d, Pln-i ^ Lp{p) > 

0. In [CF7j we showed that /3 admits a representing measure if and only if admits 
a positive extension L : V2d+2 — > Thus, as an alternative to constructing all of the 
extensions Md+i, • • • ,M2d+i, in principle it would suffice to test the positivity of the Riesz 
functional corresponding to Md+i- Unfortunately, at present there is no known concrete 
test for positivity of Riesz functionals (except in special cases, cf. |CF7| . [FNlj . |FN2j ). so 
the moment matrix extension algorithm remains the most viable approach to resolving the 
existence of a representing measure in the bivariate RD case. 

For the proof of Theorem 12.51 we require some preliminaries. For d > 1, suppose 
xi, ■ ■ ■ ,Xd are distinct and yi, . . . ,yd are distinct. Let P{x, y) := {x — xi) ■ ■ ■ {x — Xd), 
Q{x,y) := (y-yi) ■■■ {y-yd), and set Zp^Q := {{xi,yj)}i<ij<d, the common zeros of P and Q. 
Let J be an ideal in W[x, y] with real variety V = V(J) := {{x, y) S : s{x, y) = Vs G J}. 
Let /(V) = {/ G R[x, y] : f\V = 0}. In general, /(V( J)) may be strictly larger than J [CLO] . 
However, for J := {P,Q) (with P and Q as above), we will show below (Proposition 12.11]) 
that each element of /(V( J)) admits a "degree-bounded" representation which displays it as a 
member of J; in particular, J is a real ideal in the sense of [M]- Although this result may well 
be known, we could not find a reference, so we include a proof for the sake of completeness. 
First, we need three auxiliary results. 

Lemma 2.7. (The Division Algorithm in M[xi, • • • ,Xn] \CLO[ Section 2.3, Theorem 3]) Fix 
a monomial order > on Z"q and let F = (/i, • • • , fg) he an ordered s-tuple of polynomials in 
M[xi, • • • , Xn] ■ Then every f G • • • , a;„] can he written as 

f = ai/i H \- a-sfs + r, 

where Oj, G M[xi, • • • ,Xn], and either r = or r is a linear combination, with coefficients in 
M, of monomials, none of which is divisible hy any of the leading terms in fi, - ■ ■ , fg. 
Furthermore, if Oifi ^ 0, then we have multideg (/) > multideg (oifi). 

Lemma 2.8. ^ p. 67] For N > 1 let vi, ■ ■ ■ ,vn be distinct points in E?, and consider 
the multivariable Vandermonde matrix Vn '■= {vf)i<i<N,a&z'^^,\a\<N-i> of size N x .^XL^^±ii_ 
Then the rank of equals N. 
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Corollary 2.9. Let x = {xi, . . . , Xm} and y = {y\, . . . , y„} 6e sets of distinct real numbers, 
and consider the grid x x y := {(xj, yj)}i<j<m,i<j<n consisting of N := mn distinct points 
in M?. Then the generalized Vandermonde matrix Vaxy, obtained from VXr by removing all 
columns indexed by monomials divisible by or y", is invertible. 

Proof. The columns of Vn are indexed by the monomials in x and y of degree at most N, 
listed in degree-lexicographic order. The size of is N x Mi^+ii^ and by Lemma 12.81 
we know that its rank is N. We will show that T4;xy has exactly N columns, and that 
each column that was removed from Vn to produce Vx_xy is a linear combination of other 
columns in V/v. Toward the first assertion, assume without loss of generality that m < n, 
let k := n — m (so that m + k = n), and observe that the columns of V^xy are indexed by 
the following monomials: 



x,y 
x^,xy,t/^ 

. . . -r-*""! . . . 

X— ly,... ,xy""i,y™ 

m-l 2 _ _ _ m „,m+l 
y -I 1 -^y ) y 

m-l 3 _ _ _ ^„,m+l m+2 
y ) ) -^i/ 5 y 



m-l fc _ _ _ „„,m+fc-2 m+fc-1 

i/ 5 5 -'^y ) y 

m-l fc+l _ _ _ ^„,m+A;-2 
X y 1 , xy 



x™-^y""\ 

The number of monomials is then (1 + 2 + • • • + m) + mk + [(m — 1) + (m — 2) + • • • + 2 + 1] = 

'"^^"'"'^^ + m/c + ('"~^)'^ = 777,2 _j_ _ 777^j77 _|_ _ follows that Vy_xy has exactly 

N = mn columns. 

To prove the second assertion, observe that the polynomials P : — (x — xi) • • • (x — Xm) 
and Q := {y — yi) ■ ■ ■ {y — yn) vanish identically on x x y, and therefore the columns of 
Vn indexed by multiples of x"^ or are linear combinations of columns preceding them in 
degree-lexicographic order. 

By combining the preceding two assertions, it follows that Vy^xy^ having size N and rank 
A^, must be invertible. □ 



The following result is a special case of Alon's Combinatorial Nullstellensatz 
completeness, we give a proof based on Corollarv 12.91 



1; for 



Corollary 2.10. Let G = xxy be a grid as in Corollarv \2.9[ let N := mn, and letp G IR[x, y] 
be such that deg^ p < m and degy p < n. Assume also that p\g = 0. Then p = 0. 



Proof. We wish to apply Corollarv 12.91 From the hypotheses, it is straightforward to verify 
that p does not contain any monomials divisible by x™ or y", so p, properly extended with 
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zeros to indicate the absence of relevant monomials, can be regarded as a vector in M'^, the 
domain of the generalized Vandermonde matrix Vg in Corollary 12.91 Since, by assumption, 
p{xi, Hj) = for all 1 < i < m and 1 < j < n, it follows that Vcp = 0. Since Vq is invertible 
(by Corollary 12. 9p . we must have p = 0, so p = 0, as desired. □ 

Proposition 2.11. Let P{x,y) := (x — xi) ■ ■ ■ (x — Xd) and let Q{x,y) := (y — yi) ■ ■ ■ (y — yd)- 
If p := multideg (/) > d and /|V((P, Q)) = 0, then there exists u,v £ Vp-d such that 
f = uP + vQ. 

Proof. Let V := V{{P,Q)). By Lemma 12.71 we can write / = uP + vQ + r, where 
multideg (uP) < p and multideg (vQ) < p. It follows that u,v G Vp-d and that r|V = 0. 
Moreover, r is a linear combination, with coefficients in M, of monomials, none of which is 
divisible by any of the leading terms in P and Q, that is, they are not divisible by x'^ and 
y'^. Therefore, r satisfies the hypotheses of Corollary 12.101 with m = n = d. By Corollary 
[ZTnl r = 0. Thus, f = uP + vQ, as desired. □ 

Proof of Theorem \2.5[ At several points of the proof we will use the fact that if a moment 
matrix admits a representing measure and f £ Vk, then /|supp = if and only 
if f{X,Y) = in \CF2\ Proposition 3.1]. Let xi,...,Xd and yi,...,yd be sets of 
distinct real numbers, and let G := x x y = {xi,yj)^^ - denote the corresponding grid. 
Let p denote a measure whose support is precisely equal to G and let Md := Md[p\. Let 
P(x, y) := {x - xi) ■ ■ ■ {x - Xd) and let Q{x, y) := (y - yi) • • • (y - yd)- Since P\g = and 
Q\g = 0, then P{X,Y) = and Q{X,Y) = in Ca/,, whence X'^ = p{X) and = q{Y) 
for certain p,q £ Vd-i satisfying P{x,y) = x'^ — p{x) and Q{x,y) = y'^ — q{y); thus, Md is 
recursively determinate. We first show that the only column dependence relations in Md 
arise from the above relations via linearity, so that Md falls within the scope of Theorem 12.31 
If deg f = d and f{X,Y) = in Col Md, then f\G = 0, so Proposition 12.111 implies that 
there exists scalars u and v such that / = uP + vQ. Thus, f{X, Y) = uP{X, Y) + vQ{X, Y). 
Further, if deg f < d and f{X,Y) = 0, then since /|g = 0, it follows from Corollary 12.101 
that / = (whence Md-i >- 0). Thus, Md satisfies the conditions of Theorem 12.31 

Since Md has the finitely atomic representing measure p, Md admits successive positive, 
recursively generated extensions Md+i[p], Md+2[lA^ . . ., so clearly these are the unique succes- 
sive positive, recursively determined extensions of Md] let Md+k '■= Md-\-k[p] {I < k < d — 1). 
We seek to show that each of Md+i, ■ ■ ■ , M2d-i falls within the scope of Theorem 12.31 and 
that the first flat extension in this sequence occurs with rank M2d-i = rank M2d-2- We first 
give a concrete description of ker Md+k- Since Md-i 0, if r G Vd+k with f G ker Md+k-, 
then deg r = d + j ioi d < j < k. Since /i is a representing measure for Md+k, it follows that 
r I supp p = 0. Proposition 12. 1 11 now implies that there exist u,v £ Vj such that r = uP + vQ 
(with P and Q defined above in the description of p). Thus ker Md+k is indexed by the recur- 
sively determined columns; precisely, ker Md+k is the span of all of the columns x^y''{x'^ — p) 
and x^y^{y^ — q) {s,t > 0, s + t < k). Thus, Md+k satisfies the conditions of Theorem 12.31 
In passing from Md+k-i to Md+k there are d + k + 1 new columns, of which 2[k + 1) are 
recursively determined, and since these correspond (as just above) to elements of ker Md+k, 
we have rank Md+k = rank Md+k-i + {d + k + 1) — 2{k -|- 1) = rank Md+k-i + d — k — 1. 
Thus the first flat extension occurs when k = d — 1, in passing from M2d'~2 to M2d~i- D 
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We continue with an example which shows that Theorem 12.11 is no longer valid if we 
permit column dependence relations in in addition to those in (12. 3p - (12. 4p . 



Example 2.12. We define M3 by setting /3oo = ho = 
P12 = /540 = 2; /531 = /3i3 = 0; ^2 = 5, /3o4 = 22; /?5o = - 

/3o6 



= 1; Ai 

1, /341 = 



■Hm = 

-2, /332 



/3i4 = /^05 = Pm = 178, /351 = 139, (5^2 = 159, ^533 = ^W-, /324 



7 = 443272376768-2742712830r-4826809r^ _ ^hus, We have 



/?21 — /5o3 
= 13, /323 
15 



4298 



0; 
3, 



M3 



/I 








1 





1 











1 











2 


2 











1 





2 








5 


1 








2 





5 


-1 


-2 








2 





5 





-2 


13 


1 


2 





5 





22 


13 


3 





2 





-1 


-2 


13 


178 


139 








5 


-2 


13 


3 


139 


159 


2 


5 





13 


3 


894 

13 
336 

13 


159 


1657 

13 
4298 

13 


\o 





22 


3 


894 
13 


1657 
13 



2 
5 


13 

3 

894 
13 

159 

1657 

42^8 
13 
r 






\ 







22 




3 




894 








1^^7 








42^8 




13 




r 




7 


/ 



(2.10) 



It is straightforward to check that M3 is positive, recursively generated, and recursively 
determinate, with M2 >- 0, rank M3 = 7 and column dependence relations 



= p{X, y) := 40 • 1 - 2AX + 4Y - 53^^ - 2XY + 13^^ 
X'^Y = t{X, y) := 35 • 1 - 22X -Y - A6X^ + 3XY + liy^ 



and 



y^ = q{X,Y) := d 



l + d2X + dsY + diX^ + d^XY + dfsY^ + djXY'^, 



(2.11) 
(2.12) 

(2.13) 



where di 



3(-342075+1157r) 



3(487658-1651r) , 

1447 ' '^^2 — 1447 

J _ 2338519-6591r ^ _ 2(--316575+1079r) 
"5 — laaii 1 "e 



d^ 



2(-2131598+6591r) 



d7 



-48015- 



18811 

-169^ Thus, 



dA 



-2000094+6773r 
1447 ' 



M3 satisfies all of the 



18811 ' ""^ 1447 ' 1447 

hypotheses of Theorem 12. 11 except that (|2.12p is an "extra" dependence relation (not a linear 
combination of the relations defined in (j2.1ip and (j2.13p . We claim that M3 does not admit a 
moment matrix extension block -6(4) such that (M3 -6(4)) is recursively generated. Indeed, 
if such a block existed, then in the column space of (-M3 -6(4)) we would have X^Y = 
{yp){XY) := 40y - 24XY + 4y2 - bSX'^Y - 2XY'^ + 13Y^ and also X^y = {xt){X, Y) := 
35X - 22X2 _ XY - 46X3 + 3X2y + lixy^. A calculation shows that {{yp){X,Y) - 
(xt)(X,y),Xy2) = -49462+i69r ^ ^ _^ 49462^ ^3y well-defined. Thus, the 

conclusions of Theorem 12.11 do not hold for M3 (and thus there is no representing measure) . 

□ 

By contrast with the preceding example, we next show that if S RD, with all column 
dependence relations of strictly lower degree, then does admit an RG extension. 

Theorem 2.13. Suppose is positive and recursively generated, and satisfies ^2.3\) - (2m )- 
If each column relation in can he expressed as X^Y^ = r{X,Y) with deg r < i + j, then 
Md admits a unique RG extension. 
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We present the proof of Theorem 12.131 in Section [5j Finally, we note that in applying 
the algorithm, Theorem 12.31 or Theorem 12.131 may apply at some extension steps, but not 
at others. Consider |FH Example 4.15], which concerns a recursively determinate M5 with 
n = m = d = 5, deg p = 5, deg q = 4. The moment matrix M5 satisfies the hypotheses 
of Theorem 12.31 (with the roles of p and q reversed). The RG extension Mg is positive 
semidefinite and satisfies the hypotheses of Theorem 12.31 The RG extension M7 is also 
positive semidefinite, but has a new column relation, X^Y^ = r{X,Y) {deg r = 6), that is 
not recursively determined from = p{X,Y) or Y^ = q{X,Y). Thus, Theorem 12.31 does 
not apply to Afy, nor does Theorem 12.131 (since deg p = 5 = n). Nevertheless, in this case, 
when the algorithm is applied to Afy, a flat extension Mg (and a measure) results. 

3 An extension sequence that fails at the second stage 

Recall that in the most important case of recursive determinacy, a positive, flat admits 
unique positive, recursively generated extensions of all orders, M^+i, . . . M^+fc, . . ., leading 
to a unique representing measure. Further, in all of the examples of |CF3j . [CFMj and 
|F1| . when a positive, recursively generated, recursively determinate fails to have a rep- 
resenting measure, it is because it fails to admit a positive, recursively generated extension 
Md+i- These results suggest the question as to whether a positive, recursively generated, 
recursively determinate Md which admits a positive, recursively generated M^+i necessarily 
admits positive, recursively generated extensions of all orders (and thus a representing mea- 
sure) [Fll Question 4.19]. In this section we provide a negative answer to this question. In 
the sequel we construct a positive, recursively generated, recursively determinate Mi{j3^^'^) 
which admits a positive, recursively generated extension M5, but such that M5 fails to admit 
a positive, recursively generated extension Mg. It then follows from the Bayer- Teichmann 
Theorem that /^^^^ has no representing measure. 

We define M4 by defining its component blocks in the decomposition 



We begin by setting /3oo = 1^20 = /3o2 = /322 = 1, l^io = Poi = Pi2 = /324 = 2, /Jgo = /?06 = 5, 
and all other moments up to degree 6 set to 0, so that 




(3.1) 



/I 



1 

1 
1 
2 

1 

2 

1 

1 
2 



10 

2 

1 

10 

10 

2 

5 

2 

2 

2 




1 




2 

2 




M3 






1 



1 








2 




2 




(3.2) 



Vo 



5/ 



We next set 



14 



5(4) 



/2 1 





5 2 

2 
a 

6 


VO 5 



2 
2 




2\ 




2 



5 





a 

hj 



(3.3) 



where f3^Q = a, /3qi = b, /3i6 = g, /3o7 = h, and all other degree 7 moments equal 0. Let 

p{x, y) := ax^ + bx'^y + 3x^ — by — 2ax — 1 



and 



y) ■■= gxy^ + hy^ + 3y^ - 2hy - gx - 1, 
so that in the column space of (M3 i3(4)), we have the relations 



X^=p{X,Y) 



and 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



y^ = g(x,y), 

and rank (M3 5(4)) = 13. 

We complete the definition of a recursively determinate M4 by extending the relations 
(13:6]) and ([^ to the columns of (5(4)^ C(4)) , leading to 



\ 



a2 + 62 


a6 


5 





4 \ 


ab 


5 





4 





5 





4 





5 





4 





5 


gh 


4 





5 


gh 


13 + g^ + hy 



(3.8) 



Since M3 ;^ (positive and invertible), we see that M4 ^ with rank 13 if and only if 
A (4) := C(4) - B(4)^M3"^5(4) ^0. In view of ([Ml) and ([321), this is equivalent to the 
positivity of the compression of A(4) to rows and columns indexed by X^Y, X'^Y'^, XY^, 
i.e., 

'1 - 52 



[A(4)]|x3 



1 









y 0. 



(3.9) 



1 



g 



Thus, if b and g satisfy 1 — 6^ > and 1 — g^ > 0, then M4 is positive, recursively generated, 
and recursively determinate, with rank M4 = 13, so M4 satisfies the hypotheses of Theorem 

We next seek to extend M4 to a positive and recursively generated M5. In view of (13. 6p 
and (|3.7p . this can only be accomplished by defining 



X' 



{xp){X,Y) 



(3.10) 



15 



and 

■.= {yq){X,Y). (3.11) 

Theorem 12.11 implies that the resulting B{b) is well-defined and satisfies Ran -6(5) C Ran M4, 
so there exists W satisfying 5(5) = M4W. A calculation now shows that if we define C(5) 
via (I3.10p and (j3.1ip (as we must to preserve recursiveness), then M5 ^ if and only if 



A(5) = C(5) - B{5fW 



Thus, using nested determinants, and since 6^ 
generated, with rank M5 = 15 if and only if 

















o\ 


























-1+262 

-1 + 62 


bg 














bg 


-l+2g2 








-l+g2 




















Vo 














0) 



y 0. 



< 1, we see that M5 is positive and recursively 



b'< 



(3.12) 



and 



1-26^ 



2g^ + Sb^g' + bV + b'g' - b'g' 

For example, setting b = g = \, the expression in (j3.13p equals 55535 
that M5 is positive and recursively generated, with rank M5 = 15, whence M5 satisfies the 
conditions of Theorem 12.11 



> 0. (3.13) 
^995n> 0), so it follows 



With these values for h and g (or using other appropriate values), we next attempt to 
define a positive and recursively generated extension Mg. This can only be done by defining 

:= {x^p){X,Y) and := {y'^q){X,Y). Theorem O implies that the resulting B{6) is 
well-defined and that there is a matrix V such that -6(6) = M^V. Further, C(6) is uniquely 
defined via the preceding column relations. Mq as thus defined is recursively generated (by 
construction), but we will show that it need not be positive. Indeed, a calculation shows 
that A(6) = C(6) — -B(6)^y is identically except perhaps for the element in the row and 
column indexed by X^Y^ (the row 4, column 4 element), which is equal to 

(1 - 36^ + b^- ab^g + ab^g + bh- 2b^h){-l - ag + 3g^ + 2ag^ - + bg^h - bg% 
-1 + 262 + - 362^2 _ 54^2 _ 52^4 + 54^4 • 

Note that the denominator of the preceding expression is the negative of the expression in 
(j3.13p . and is thus strictly negative. Thus Mq is positive if and only if 

7] := {l-2,b^+b^-ab^g+ab^g+bh-2b^h){-l-ag+'ig^ = 2ag'^ -g'^+bg^h-bg^h) < 0. (3.14) 

With 5 = (7 = i, we have 

_ (-836 + 15a - 224/1) (836 + 224a - 15/i) 
^ ~ 1048576 ■ 

If we choose a and h so that 77 = 0, then Afg is a fiat extension of M^, and /3 = /S^^^ has 
a 15-atomic representing measure. If we choose a and h so that < 0, then Mq is positive 
with rank 16, and since, in Corollary 12.41 = m = 4 and d = 6, it follows that Mq has a fiat 



extension M7. However, if we choose a and h so that 77 > (e.g., with h = and a > ^), 
then Mq is not positive, whence /3 has no representing measure. 
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4 Proof of Theorem 12.11 



The proof of Theorem 12.11 entails two mains steps: (i) the construction of the block B{d + 1) 
from the column relations ()2.ip and (j2.2p so that (Md B{d + 1)) is recursively generated; and 
(ii) the verification that Ran B{d+\) C Ran M^. 

STEP (i): Step (i) will follow from a series of five auxiliary results (Lemmas 14 . 1 1 - . 
To begin the formal definition oi B[d+\), note that blocks i?[0, d+ 1], . . . , -B[(i — 1, d+ 1] are 
completely defined in terms of moments in M^. Indeed, for 0<i<(i+l, 0<j<(i — 1, and 
h,k > with h + k = j , the component of B[j, d + 1] in row X'^Y'' and column x^y'^+^-\ 
which we denote by {X^Y'^^^~^ , X'^Y'') , must equal Pi^h,d+i-i+k- Note also that for i > n, 
the above component is alternately defined by (I2.8p . so we must show that the two definitions 
agree. 

Lemma 4.1. For 0<f<d+l — n and i,j > with i + j < d — 1, the entry in column 
^n+fyd+i-n-f^ j-^y; x^Y^ , as defined by I12.8\) . coincides with the moment inherited from 
Mfi by moment matrix structure, j5n+f+i4-n-f+j+\- 
Proof. Consider first the case when d — n — f > 0. From (j2.8p . we have 

= Yl arsX^'+^Y'+'^+^-'^-f (0 < / < d + 1 - n), 

r,s>0,r+s<n—l 

SO 

Since r + f + s + d+1 — n — f<d, s + d+ 1 — n — f > 1 and i + j < d—1, using the moment 
matrix structure of the blocks of M^i we may express the last sum as 

Ors (X'^+Zy «+'^-"-/ ,X'Y^+^). 

Now ()2.3p implies that in the later expression is equal to 

(x"+/y'^-"-/,x*y^-+^) = 

For the case f = d — n + 1 and i + j < d — 1, 

{X'^+\X'Y^) = 'Y^rs{X''Y'X'^+^-'',X'Y^) 
= Y^rs{X''Y'X'^~'',X'+^Y^) 
= {X'^,X'+'^Y^) 

= /3d+i+lJ- 

□ 

We have just verified that in the left recursive band, in blocks of degree at most d—1, 
each column element coincides with the corresponding "old" moment from M^. Old moments 
are also used to define the central (nonrecursive) band of columns in blocks of degree at most 
d—1. We next use these left and central bands, together with ()2.9p . to show that the 
column elements in the right recursive band, in blocks of degree at most d — 1, also agree 
with corresponding old moments. 
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Lemma 4.2. For 0<k<d + l — m, i,j>0,i + j<d— 1, column X^'^^ kym+k ^ 
defined by satisfies {X'^+i-m-kym+k ^ x'Y^) = (3,+d+i~^-k,m+k+j ■ 

Proof. The proof is by induction on k. For A; = 0, we show that (^X'^~^^~'^Y"^, X^Y^) = 
Pi+d+i-m,m+j ■ From ([IJ]), we have 

= Yl aab{X'^^^-"'+''Y\X'Y^) 

a,b>0,a+b<m-~l 

+ Yl -/ce{X^^'^+'^~"'Y^,X'Y^). 

c,e>0,c+e=m,e<m 

Since a + b < m, then in M^, 

Since e < m, (^X'^^'^^^~'^Y^ , X'^Y^) is in either the left or central band, and thus equals the 
old moment /3c+d+i-m+i,e+i- Now 

—'m+a+i,b+j + ^ ^ Tce/^c+d+l— m+i,e+j • 

a,b>0,a+b<m—l c,e>0,c+e=m,e<m 

In Md, the latter expression equals 

as desired. We next assume the result is true for 0, . . . ,k — 1. Consider first the case 
when k < d + 1 — m. We have 

^Xd+i-m-kym+k^X^Y^) = J2 aab{X'^^^-'^-''+''Y''+\X'Y^) 

a,b>0,a+b<7n~l 

+ Y ^ce{X^^'^+^~'^-''Y''+'',X'Y^). 

c,e>0,c+e=m,e<m 

The term (^X'^~^^~"^~^^"'Y^~^'' , X'^Y^) is a component of M^, and thus equals the corresponding 
moment. Since e + k<m + {k — 1), X'^'^'^~^^~"^~^Y^^^ is, by induction, a column for which 
the elements of row-degree i + j are old moments. Thus, 

^Xd+l-m-kym+k ^ X^Y^) = OiabPd+l-m-k+a+i,b+k+j + 7cePc+d+l-m-k+i,e+k+j ■ 

In Md, the last expression equals 

_|_ ^ ^ ^j^c+d-m-fc-y^e+fc X^^^Y-^) = (^X'^~'^~^Y^~^^ X^~^^Y-^) 



Pd—m—k+i+l,m+k 



+3- 
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Finally, we consider the case k = d + 1 — m. We have 

a,b>0,a+b<m—l 

+ -fce{X''Y^+'^+^-"',X'Y^). 

c,e>0,c+e=m,e<m 

Since e < m, then c > 1, so X'^Y'^~^'^~^^~"^ is to the left of i.e., c = d + 1 — m — k' for 

k' = d+ l— m — c<k. Thus, by induction, 

{Y''+\X'Y^)=Y^ab|3a+^,b j + 7ce/3c+i,e+d+l-m+i- 

In Mrf, the last expression equals 

J2aab{X^Y^+'^-'^,X'Y^+^) +Y^ce{X''Y^+'^-'^,X'Y^+^) = {Y'^,X'Y^+^) 

= (3i,d+j+l, 

as desired. □ 

To complete the definition of B(d + 1) we must define B[d, d + 1]. Within this proposed 
block, we first use (|2.8p to define the left recursive band, . . . , x^Y'^~^^~"' . Note that be- 

tween the end of the left band, x^Y^~^^~'^ , and the beginning of the right band, x"-^^~"''-Y"^ , 
there is a central band oin + m — d — 2 columns; set 6 := n + m — d — 1. In row X'^, 
each of the components in the central columns, (x"-iy'^+2-", X^), . . . , (x'^+^-^y'""!, X^), 
corresponds via a cross-diagonal to a component of column X"y^+^~" (whose value is known 
from dlSD), i.e.. 

We may thus use (j2.9p to define X'^), and we extend the latter value along 

the central-band section of the cross-diagonal to which it belongs. Next, in row X'^~^Y, 
we use this value with (j2.9p to define {X'^~^^~"^Y"^ , X'^~'^Y) , and we extend this value along 
the central-band section of its cross-diagonal. Proceeding in this way, we completely define 
column X'^+i-'"y"" and insure that it is Hankel with respect to the central band. Finally, 
we use dZl]) to define column X'^-^F'^+i, and, successively, X'^-"^-'^Y"'+^ , Y'^+^. This 
completes the definition of a proposed block B[d,d+ 1]. However, to ensure that it is well- 
defined as a moment block, we must check that for a cross-diagonal which intersects columns 
^nyd+i-n ^^^^ j^d+i-mym^ ^^le Components of the cross-diagonal in these columns agree in 
value, i.e., the values arising from (j2.8p are consistent with those arising from ()2.9p . More 
generally, we need to show that the block we have defined is constant on cross-diagonals. 

To show that B[d,d + 1] is well-defined and Hankel, we begin with the following general 
result concerning adjacent columns that are recursively determined from the same column 
dependence relation. Suppose in Col there is a dependence relation X^Y*^ = p{X,Y)^ 

where c + e = d and p{x,y) = ^2 '^abx"'y^ G Vd-i- Then the elements of Col 

a,b>0,a+b<d-l 

defined by 

X'+'Y' ^ {xp){X,Y) := Yl c^abX^'+'Y' 

a,b>0,a+b<d-l 
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and 



X^y^+i = iyp){X, Y) := Yl c.abX''Y'+' 

a,b>0,a+b<d-l 

are Hankel with respect to each other, as fohows. 
Lemma 4.3. For i,j > 0, i + j < d, j > 0, 

(X'^+^Y^X'Y^) = {X^Y^+\ X'+^Y^-'^) 

Proof. We have 

(^X^+^Y^^X'Y^) = Yl aabiX'^+^Y^X'Y^), 

a,b>0,a+b<d-l 

and since each row and column in the last sum has degree at most d, relative to we may 
rewrite this sum as 

Y aab{X''Y''+\X'+^Y^-^) = {X^Y^+\X'+^Y^-^). 

a,b>0,a+b<d-l 

This completes the proof. □ 

It follows immediately from Lemma [4.3l that the left recursive band in B[d, d+1] is constant 
on cross-diagonals. We next check that if an element of a column in the non-recursive central 
band can be reached on a cross-diagonal which intersects both columns X"^Y'^^^~'^ (at the 
edge of the left recursive band) and X'^~^^~^Y^ (at the edge of the right recursive band), 
then the values obtained from both of these columns agree. This is the substance of the 
following lemma. 

Lemma 4.4. For Q < k < 2d + 1 — n — m, 

Proof. The proof is by induction on k. We begin with the base case, A; = 0, and seek to show 
that = (x"y'^+i-",X'^-''y'5) (recaU that 6 := n + m - d - 1). Using 

(IZQl) . we may express (X'^+i^^y'", X'^) as 

Y aab{X'^^^-'^+''Y\ X'^) + Y ice{X'^^^-^Y'',X'^). (4.2) 

a,b>0,a+b<m—l c,e>0,c+e=m,e<m 

Note that (^x'^+i-m+ayb ^ x'^) is a component of M^; further, since e < m, (^x^+'^+^-'^Y'' , X'^) 
is the endpoint of a cross-diagonal that lies entirely in the left and central bands, and is thus 
constant. Therefore, we may rewrite (|4.2p as 



a,b>0,a+b<m—l c,e>0,c+e=m,e<m 
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= ars Y aafe(^'^""+'T^ X'^+l-^+'^y'') + Y "rs J]] 7ce (^"'""+''+^5^" , ^'^"^1^') 

Since 5 = m + n — d— 1, in the last sum is equal to 

which completes the proof of the base case. 

We assume now that ()4.ip holds for 0, . . . , A; — 1, with k — l<2d+l — n — m. To establish 
()4.ip for A:, we consider first the case d — k > n. Let us write k := (^X'^~^^~'^Y"^ , X'^~''Y'') as 

K = ^ Q„fe(x'^+l-'"+'^y^x'^-'=y'=) 

a,b>0,a+b<m,~l 

+ ^ 7ce(x'^+^"^y',x'^-'=y'=) 

c,e>0,c+e=m,e<m,d+l— e>n 

+ Y 7ee(X'^+^"'y%X'^~^y''>. (4.3) 

c,e>0,c+e=m,e<m,(i+l— e>n,ci+l— e<n 

Note that the components in the first sum of (j4.3p lie in M^. In the third sum, since d+l — e < 
n, column X'^^^~^Y'^ is in the middle band, and the component 7 := (^X'^~^^~^Y'^ , X'^~''Y^) 
lies on a cross-diagonal a strictly above the cross-diagonal for k. Either because a does not 
intersect column x'^+i-™y™^ or by induction if it does, we see that 7 has the same value as 
^j^nyd+i-n^^d-fc-(n-(d+i-e))yfc+n-(d+i-e)^ ^j^g ^^^^ cross-diagonal) . Thus gS]) Can 
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be expressed as 

c,e>0,c+e=m,e<m,ii+l— e>n 
_j_ ^ ^ ^ ^j^n^d+l-n j^d-fc-{n-(d+l-e))-y^fc+n— e)^ 

c,e>0,c+e=m,e<m,d+l— e>n,d+l— e<n 
_|_ ^ ^ ^ ^j^n-^d+l— n j^d— fe— (n— e))-y^fc+n— (d+l— e)^ 

= ^ Oaf, 5^ a^,(x^+'^-^-"y^+^ 

_|_ ^ ^ ^ ^j^rys+d+l-n j^d-fc-(n-(d+l-e))yfc+n-((i+l-e)^-j (4 4) 

Using the symmetry of in the first and third inner sums of the last expression, we may 
rewrite this expression as 

_|_ ^ ^ ^ ^j^a!-fc-(n-(d+l-e))yfc+n-(d+l-e) j^rys+d+l-n^ -j (4 5) 

In the second inner sum of ()4.5p . (^X'^~^'^~^^~^~"'Y'^~^^ , X'^~''Y^) is a component of and 
thus equals the moment I3r+d+i-e-n+d-k,s+e+k- Since X'^~''~^^~'^Y'^~^^ is a row of degree at 
most d—1, this moment coincides with (^x'^~^'^~^^~"^Y'^ , X'^~''~^'^""Y^~^'') from the left band of 
B[d + r + s — n,n + 1]. Further, in the third inner sum of ()4.5p . 



is also a component of Md, equal to /3r+d-k-{n-{d+i-e)),s+k+n-{d+i-e)+d+i-n, and this mo- 
ment coincides with {X'^~^^~'^Y'^ , x'"+'^~'=~"y*+'=) from the middle band in B[d+r+s-n, d+1]. 
Thus, the expression in (j4.5|) can be written as 

Y^^rsi Yl aab{X'^^^~'^+''Y\X''+'^-''-''Y'+'') 

a,b>0,a+b<m—l 

_j_ ^ ^ ^ m+c-^e j^r+d— fc— n-^s+A.'^ 

c,e>0,c+e='m,e<m,,d+l—e>n 

+ Y Y 7ce(x'^+^"™+'y',x^+'^-^-"y''+*^)), (4.6) 

c,e>0,c+e=m,e<m,d+l—e<n 
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which equals 

a„(X'^+i-™y'", x^+<i-fc-«y^+fc)). (4.7) 

Since x^~^'^~^~"'Y'^~^'^ is a row of degree at most d — 1, Lemma [4.21 imphes that the expression 
in (|4.7p equals 

^ ^ 0'rsPd+l—m.+r+d—k—n,m.+s+k — ^ ^ '^rs/^r+d— (5— fc,s+a!+l— n+(5+fc 

This completes the proof of the induction step for (j4.ip when d — k > n. 

We next treat the case when d — k < n, which implies 5 + k > m. We have 

a,6>0,a+t'<m— 1 
c,e>0,c+e=m,e<m 

(4.8) 



Note for future reference that all of the matrix components that appear in (j4.8p come from 
Md. 

We now consider 

a,b>0,a+fe<m— 1 

c,e>0,c+e=m,e<m 

= ^a„b(x'^-'^y^x'^+l-™+"y*) 

+ ^7ce(x'^+^"'y^x'^"''y'') (4.9) 

(using symmetry of Al^ in the first sum). Since k — {n — (d — k)), d+1 — m + a — n + d— k 
(= a + d — d — k), and h + n — {d — k) are all nonnegative, by applying the block-Hankel 
property of M^, to the first sum in (j4.9p . we may rewrite the expression in (j4.9p as 

+ ^7ee(x'^+^"'y^x^"''y^') (4.10) 

^ ^ ^ ^ ^ ^ ^ ^j^ryAS+fc— (n— (d— fe)) j^d+1— -m+a— n+d— (d— fc) ^ 

+ J]]7ce(x'^+l-^y^x'^-'=y'=), (4.11) 
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and all of the matrix components in the first double smn of (j4.1ip are from M^. Comparing 
the components in the first double sums of ()4.8|) and ()4.1ip . we have 

/ \ra+d~5—k\/-b+5+k-~in \rr\^d+l~n+s\ o 

\^ -I ) ^ J / — Pa+d-5-k+r,b+5+k-m+d+l-n+s 

= l^r+d+l-m+a-n+d-k,s+k-{n-(d-k))+b+n-{d-k) 

^j^r-^s+fc— (n— (d— fc)) m+a— n+d— fc-^b+n— (d— fc) ^ 

SO the first double sums of (14.81) and ()4.1ip are equal. 
Let us write the rightmost sum in (j4.1ip as 

c,e>0,c+e=m,e<m,d+l—e>n 

+ Yl -fce{X'^^^-^Y^,X'^-''Y''). (4.12) 

c,e>0,c+e=m,e<m,d+l—e<n 

In the second sum of (j4.12p , since d+l-e < n, the component (^x'^+i-^y'' , X'^-'^V'') (from the 
middle band) has the same value as the component (x^y^+i"", X'^-'=^("-('^+^+^)y'^'+"^('^+i^^)) 
on the same cross-diagonal. (This is because the cross-diagonal is strictly above that for 
^^d+i-mym^ j^(i-fcyfc^^ SO the conclusion follows by definition or induction.) We may now 
write the expression in (j4.12p as 

c,e>0,c+e=m,e<m,d+l—e>n 
I ^ ^ ^ j^nyd+l— n j^d— fc— (n— e))-yA;+7i— (d+l— e)^ 

c,e>0,e<m,c+e=m,d+l— e<n 
c,e>0,e<m,c+e=m,d+l— e>n 

+ Y 7ceXl«rs(^''^'+'^+^"'',^''"*'"^""^'^+^"'^^l'''+""^'^+^"'^)-(4.13) 

c,e>0,c+e=m,d+l—e<n 

All of the matrix components in (j4.13p are from M^, so (j4.13p can be expressed as 

c+e=m 

It is straightforward to check that this double sum coincides with the second double sum in 
()4.8p (whose matrix components also come entirely from M^). This completes the proof that 
the second double sums in (j4.8p and ()4.1ip have the same value, so the expressions in (j4.8p 
and (14. lip are equal, which completes the proof of the induction when d — k < n. Thus, the 
induction is complete. □ 

We have shown above that in B[d,d + 1] the columns X'^'^^, . . . ^ x'^~^^~"^Y"^ are well- 
defined and Hankel with respect to one another. Using ()2.9p . we also successively defined 
columns . . . , Y'^+'^. We next show that the columns X'^-'^+^Y"', F^+i are 

Hankel with respect to each other, so that all of B[d,d+ 1] has the Hankel property. 
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Lemma 4.5. For 0<s<d + l — m and i,j > with i + j = d and j > 0, we have 

Proof. The proof is by induction on s > 0. For s = 0, we have (^X'^~^^~"^Y'^ , X^Y^) 

J2 aafe(X'^+^"™+''y^X*y^■) + -fce{X'^^'^-^Y^,X'Y^). (4.14) 

a,b>0,a+b<m~l c,e>0,c+e=m,e<m 

In the first sum, each component is from M^. In the second sum, column X'^~^^~'^Y^ is strictly 
to the left of X'^~^^~'^Y"^ , so it is Hankel with respect to its right successor, 
may thus rewrite the expression in (j4.14p as 

aab{X'^~"'+''Y''+\ X'+'^Y^-'^) + Y jce{X'^~^Y^^\ X'+^Y^^'^) 

a,b>0,a+b<rn—l c,e>0,c+e=m,e<m 
_ ^j^d— m^m+1 X^~^^Y-^~^) 

Assume now that the Hankel property holds through s — 1 and consider 

^^d+l-m-sym+s^j^iYj^ = ^ a^f^{X'^+'^"-"'+''~'Y''^' , X'Y^) 

a,b>0,a+b<m—l 

+ Y ^ce{X'^^^~^~'Y''+',X'Y^). (4.15) 

c,e>0,c+e=m,e<m, 

As above, in the first sum, each component is from M^; in the second sum, each column 
j^d+i-e-sye+s [g ^]^q jgf^ j^d+i-m.-sym.+s ^ qq ^j^g Hankel property holds for this column 

by induction. We may thus write the expression in M.lSh as 

Y aab{X'^~"'+''-'Y''+'+\X'+^Y^-'^) + Y -fce{X'^~^-'Y^+'+\X'+'^Y^-^) 

a,b a+b<m-l c+e=m,e<m 

which completes the proof by induction. □ 

STEP (ii): The preceding results show that under the hypotheses of Theorem 12. 11 there 
exists a unique block B{d + 1) that is consistent with recursiveness in [M^, B{d + 1)) . To 
prove Theorem 12. H we must also show that Ran B{d+ 1) C Ran M^- The following lemma 
is a step toward this end; it shows that the rows of (Mrf B{d+l)) of the form X'^+fYd 
if^g^O, n + f + g < d) are recursively determined from row X"'. 

Lemma 4.6. For i, j > 0,i + j < d + 1 and for f,g>0, n + f + g < d, 

(X'y^ X"+^yf ) = Y ars{X^Y^,X''+fY'+3), (4.16) 

r,s>0,r+s<n — 1 

Proof. Since is real symmetric, it follows from (|2.8p that (j4.16p holds for i + j < d. We 
may thus assume that j = d + 1 — i. Consider first the case when n + f -\- g < d. In the 
subcase when i < d, it follows from the presence of old moments in B\n + f + g,d + V\ that 



{X'Y^+^-\X^^fY3) = (3,+n+f,d+i- 



■i+gi 
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and in we have 

ars{X~'+fY'+3+i^X^Y'^~^) (4.17) 

r,s>0, r+s<n— 1 

= Yars{X'Y'^^\X'^+fY'+3+i) (by symmetry in Mrf) 

= Ors/5i+r+/,a!-j+s+g+l 

= Yars{X^Y'^+^-\X''+fY'+3) 

(by moment matrix structure in B[d + 1)). 

For the subcase when i = d + 1, we first note that {X'^'^^ , X^'^fY^) = f^d+i+n+f,g = 
{X'^, X"-+fY9) = {X''+fY9,X'^), and we then proceed beginning as in Uti\ . 
We next consider the case n + f + g = d, and we seek to show that 

^^iyd+i-i^ ^n+Zyc,) = ars{X'Y'^+'^-\X''+fY'+s). (4.18) 

We begin by showing that (j4.18p holds if the column x*y^+^~* is recursively determined 
from ()2.8p . i.e., i > n. In this case, we have 0<i<d+l — n, so 

(X*y'^+i"\X"+^ys) = Y,ars{X''Y'X'~''Y'^+'^-\X''+fY3) 

= ^ a„„ ^ a„(x"y''x^yf , x^-y^x^-^y^+i-*) 
= 5^ a„(x^y'^x^-"y'^+l-^ 

Thus 

{X'Y'^+'^-\X''+^Y3) = ^a™(x*y'^+i-^x"+^y''+3), 

which is equivalent to (j4.18p . 

Returning to the proof of M.lSh . we next assume that column x^Y'^~^^~'^ is not recursively 
determined, i.e., d+1 — m<i<n. By the Hankel condition in B{d +1), we have 

(^X^Y'^^^~^ X'^'^^Y^) = ^j^"-y^+i-" jf«+/y"-«+9^ 

= Z a.rsi3r+i+f,s+d+l-i+g (in Mrf) 

= J]]a,.,(xy"'+i-\x''+-^y"+'?) 

(since r + f + s + g<n + f + g = d). 

Note that if ()4.16p holds for a collection of columns, then it holds for linear combinations 
of those columns. Thus, using the preceding cases and (|2.9p . we see that (j4.16p holds, 
successively, for x'^'^^~"^Y"^ , . . . , y^"*"^, which completes the proof. □ 
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The following result shows that the rows of (M^ B{d + 1)) of the form X^I^Y'^^^ (/, g > 
0, m + f + g < d) are recursively determined from row y™. 

Lemma 4.7. For i, j > 0,i + j < d -\- 1 and for f,g>0, m + f + g < d, 

{X'Y^ ^Xfy^+a) = 6„„(xV^X"+-^y''+9). (4.19) 

u,v'>0,u+v<.m,v<m 

Proof. Since is real symmetric and recursively generated, its rows are also recursively 
generated from (j2.ip and (|2.2p . so (j4.19p holds ii i+ j < d. We may now assume j = d+l — i, 
and we first consider the case m + f + g < d and the subcase i < d. Since f + g + m < d, 
using old moments we see that 

= {X^Y'^+s+i^x'Y'^-') {in Md) 

= ^6™(x"+^y''+f+\x'y'^-*) 

= (in Md) 

= J]]6„„(x*y'^+i-\x"+-^y^+9) 

(since u + v + f + g<d). 

The subcase when i = d + 1 proceeds as above, but starting with {X'^^^ , X-I'Y™''^^) = 
Pd+i+f,m+g = {X'^,Xf+^Y"'+9) = {Xf+^Y"'+9,X'^). For the case m + f + g = d, we 
first consider the subcase when i > n, so x*y^"^^~* is in the left recursive band. We have 

= ^ars{X'-+^-''Y'+'^+^-\XfY'''+3) 

= ^arsY^ 6„^(x''+^-"y"+'^+i^\ x"+^y^+f ) 

(by row recursiveness in M^) 

= ^6™(x^y'^+l-^x"+^y''+5). 

In the next subcase, we consider a column in the center band, of the form X'^~^^~^Y^ with 
d + 1 — n<i<m. In this case, ()4.19p is equivalent to 

(^X'^+^~iY\XfY"'+3^ = buv{X'^^^^'Y\X''+fY"+3), (4.20) 

u,v>0,u+v<m,v<'m 

Note that the component (X'^+^-^y*, X^y^+s) lies on a cross-diagonal that reaches column 
j^d+i-mym^ SO since B{d+ 1) is well-defined, we have 

Y 6™(x"+'^+^-'"y^x■^+™-^y»+*). (4.21) 

u,v'>0,u+v<.m,v<m 

For the subcase when u + v < m. Md we have 

\^ r ,A-' y ) — Pu+d+l+f-i,v+g+i 

= (x'^+i-V*,x"+^y''+f) (since n + / + + 5 < d- 1). 
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For the subcase when u + v = m, there are three further subcases in showing that 
For V = i, (|4.22p is clear. For v < i, the Hankel property in B[d, d+ 1] implies 

For V > i we have, similarly, 

^j^d+l—iyi j^u+fyv+g\^ ^j^d+l~i—{v—i)yi+v—i j^u+f+v—iyv+g—{v—i)\^ 

^j^d+l+u—myv J^m+f—iyg+i^ 

Since (|4.19|) holds in and in all columns of the left and center bands, it now follows, using 
(j2.9p successively, that it holds for columns in the right recursive band, which completes the 
proof. □ 

We are now prepared to prove that Ran B{d + 1) C Ran M^. It follows immediately 
from (12. 8p that each column in the left recursive band of B{d + 1) belongs to Ran M^. In 
view of (j2.9p . to establish range inclusion, it suffices to show that each central-band column 
of B{d+ 1) belongs to Ran Md- Let S denote the set of recursively determined columns of 
Md, i.e., 

g l^^^*^ X^^^ X^Y X'^ j^nyd—n ym XY"^ ym+l j^d—mym Y"^^ 

Let B denote the basis for Col (the column space of M^) consisting of those columns of 
Md which do not belong to 5. Let Mg denote the compression of Md to the rows and columns 
indexed by B. Since Md h 0, we also have Mb h 0. Let X'y^+i"* {d + 1 - m < i < n) 
denote a central-band column of B{d+ 1), and let Vi = [X^Y'^^^~^]b denote the compression 
of x^y^+i-* to the rows of B. There exists a unique vector of coefficients {CaDx'^V'eB such 
that 

i.e., for each X"y^ G B, 

^^iyrf+l-i^^«y^^ ^ ^ c2(^"5^^^"l^'')• (4.23) 

To complete the proof that Ran B(d + 1) C Ran Md, it suffices to prove that = 
which, in view of (I4.23p . follows from the next result. 

X'^ybf^B 

Lemma 4.8. For each X'^Y^ G S, 

{X'Y'^+^^\X'^Y^) = c^^^X'^Y^X'^Y^). (4.24) 

X^Y''&B 
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Proof. We may assume without loss of generality that n < m, so the elements of S may 
be arranged in degree-lexicographic order as X", • • • , Y"^, ■ ■ ■ , X'^, ■ ■ ■ , Y'^. We will prove 
(j4.24p by induction on the position number of row X^Y'^ £ S within the degree-lexicographic 
ordering. For row X" (c = n, e = 0), Lemma 14.61 implies that 

(XV^+i^\X")= ars{X'Y'^+^~\X'"Y'). (4.25) 

r,s>0,r+s<n~l 

Since r + s < n, X'^Y^ G B, so the sum in (|4.25p may be expressed as 

(using Lemma 14.61 again) . Assume now that ()4.24p holds for all rows X^Y^ S S with order 
position up to k — 1, and consider X'^Y'^ G S with position k. Either c > n or e > m; we 
present the argument for the case e > m (the other case is simpler). We have e = m + g for 
some 5 > 0. From Lemma 14.71 we have 

u,v>0,u+v<m,v<m 

Now X'^'^^Y^'^'" is either a basis vector, or, since f < m, it precedes X'^Y^'^^ in the ordering 
of 5. Thus, by definition (for the basis rows) and by induction (for the non-basis rows), the 
preceding sum is equal to 

= E^™ E c^!{x^y',x'+^Y3+-)= c^l{x''Y\YbuvX'+^Y3+-) 

X°-Y^&B X'^Ybf^B 

= Y c^2{X'^Y\X^Y^) 

X^^Y^&B 

(by another application of Lemma 14. 7p . □ 
The proof of Theorem 12.11 is now complete. 

5 Proof of Theorem [2J^ 

For the proof of Theorem 12.131 we require a preliminary result concerning a general moment 
matrix. 

Lemma 5.1. Suppose M^+i satisfies Ran B{d+ 1) C Ran M^. If p S Vd and p{X,Y) = 
in Col Md, then p{X,Y) = in Col Md+i- 

Proof. Since is real symmetric, we have p{X, 1") = in the row space of M^, and we 
first show that p{X,Y) = holds in the row space of [Ma B{d+ 1)) . Let p := deg p and 
suppose p{x,y) = Ylr s>q r+s<p^^sX^y^ ■ Then for i,j>0 with i + j < d, we have 

Y «rs {X'Y^ ,X'-Y') = 0. (5.26) 

r,s 
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Consider a column of degree d + 1, X^y^+i-" (0 < n < fi + 1). We seek to show that 

J2ars{X''Y'^+^-'',X'Y') =0. (5.27) 

r,s 

By the range inclusion, we have a dependence relation in Col [M^ B{d + 1)) of the form 

^uyd+i-u^ c^^^X^yK (5.28) 

a,b>0,a+b<d 

Thus, 

r,s r,s a,b>0,a+b<d 

Now, p{X, y) = in the row space of (M^ B{d + 1)) , so p{X, y) = in Col (^^(^f^i^T^ ■ 

□ 

Proof of Theorem \2.13[ It follows from the proof of Theorem 12.51 that admits a unique 
extension M^+i which satisfies Ran B{d + 1) C Ran and such that (j2.8p - (|2.9p hold 
in Col Md+i- It remains only to prove that Md+i is recursively generated. Since Ma is 
recursively generated, it suffices to consider a dependence relation in Col M^+i of degree d, 
say 

^iyrf-i ^ CghX^Y^ (5.29) 

g,h>Q,g+h<d-l 

(where < i < d), and to show that 

j^i+iyd-i^ Y CghXS+^Y'' (5.30) 

g,h>Q,g+h<d~l 

and 

J^^Yd~^+l ^ Y CghX^Y^^+K (5.31) 

g,h>0,g+h<d-l 

Suppose first that i > n, so that x^Y'^~^ lies in the left band. Then from (12. 3p we also have 

XiY'^-i= Y arsX'^''+''Y'+'^-\ (5.32) 

r+s<n—l 

Thus, in we have the column relation of degree at most d — 1, 

^ CghX^Y'' = arsX'-''+''Y'+'^-\ 

g+h<d—l r+s<n—l 

Since is recursively generated, it follows that in Col we also have 

Y CghX3+^Y'^= Y arsX'-''+''+^Y'+'^-\ 

g+h<d—l r+s<n—l 
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Lemma 15.11 implies that the last equation also holds in Col M^+i, where, from (j2.9p . the 
right-hand sum represents x^~^^Y'^~^; this estabhshes (I5.30p . We omit the proof of ()5.3ip . 
which is similar. The case when d — i > m, so that x^Y'^~^^^ is in the right band, is handled 
in an entirely analogous fashion, so we also omit the proof of this case. 

We next consider the case when d — m < i < n, so that column x^Y'^~^ in ()5.29p is in the 
central band. To establish (j5.30p . it suffices to verify that 

^^i+iyd-i^^fcyj^ = Cgh{X'^+^Y'',X''Y^) {k,j > 0, k + j<d+l). (5.33) 

9,h>0,g+h<d-l 

The case when k + j < d is easy, using (j5.29p and the old moments in block B[k + j,d + 
1]. We consider next the case k + j = d and the subcase when k > n. In this sub- 
case, (^x^~^^Y'^~^ , X^Y'^~^) belongs to a cross-diagonal of B[d,d + 1] that intersects column 
^nyd+i-n^ SO from the definition of B[d,d + 1] in the proof of Theorem 12. H we have 

= ^a,.,(X^-y"+"'+i-",X'=-("-*-i)y°'-'=+"-*-i). (5.34) 

Now, we have 

^ Cgh.{x3+^Y^,x''Y'^~'') = ^c3/,(x'=y"'-^x3+ly^) 



c, 

g,h>0,g+h<d-l 



Cgh Y «rs(X''+''~"l"'+'^"^ X3+^Y^) 
rs 

Y ars Y CghiX'^^^Y'', ;^r+fc-ny.+d-fc^ 

Y E ^9h{X'^Y\ x'^+'=-"+iy^+^-'=) (in Md) 

Y ars{X'Y'^~\X''+''-''Y'+'^-'') 

Y ar•s(^'■+''""i^'+"'"^ X'Y'^-') 



d—i—{k—n+l) \ 



This last expression agrees with (j5.34p . so (j5.30p is established for this subcase. The proof 
of this subcase for (I5.3ip is very similar, so we omit the details. In the subcase when k < n, 
thend— A; > m, and we see that (^x^~^^Y'^~^ , X'^Y'^~^) belongs to a cross-diagonal of i?[(i,(i+l] 
that intersects column x'^~^^~'^Y'^ . Since deg q < m, the proof of this subcase is entirely 
analogous to that above, but using (j2.9p for the definition of X'^~^^^™ Y"^ . 

Finally, we consider the case k + j = d + 1. As above, we will treat the subcase of (j5.30p 
when k > n in detail and omit the proofs of the other subcases of (j5.30p and (|5.3ip . which 
are similar. Since k > n, then, as above, we have 



hn— i— 1\ 



= ^a^,(X^y^+'^+i-'^,X'=-("-'-i)y'^-'=+"-*). (5.35) 



Now, 
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g,h>0,g+h<d~l 

(since (M^ B{d + 1)) is the transpose of 



rs 

Since the row degrees of the terms in the last sum are at most d, by the previous cases (for 
j + k < d and j + k = d), the last double sum may be expressed as 

relative to {Ma B{d + 1)) . Since M^+i is real symmetric, the latter sum may be expressed 
as 

Yars{X''+^-'^Y'+'^+^-^,X'+^Y'^-') 

^ ^ ^ f^^r Y s+d+l—n) ^i+l+k~nYd~i—k+n\^ 

and this agrees with ()5.35p . The proof is now complete. □ 
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